The unification of inflation and late-time acceleration in the frame of /c-essence 
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By using the formulation of the reconstruction, we explicitly construct models of fc-essence, which 
unify the inflation in the early universe and the late accelerating expansion of the present universe 
by a single scalar field. Due to the higher derivative terms, the solution describing the unification 
can be stable in the space of solutions, which makes the restriction for the initial condition relaxed. 
The higher derivative terms also eliminate tachyon. Therefore we can construct a model describing 
the time development, which cannot be realized by a usual inflaton or quintessence models of the 
canonical scalar field due to the instability or the existence of tachyon. We also propose a mechanism 
of the reheating by the quantum effects coming from the variation of the energy density of the scalar 
field. 
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INTRODUCTION 



We now believe the accelerating expansion of the present universe by several cosmological observations |lHj]. The 
acceleration has been often supposed to be generated by the dark energy, which is an unknown fluid. So-called k- 
essence model 0-01 is a model of the dark energy. The fc-essence model is originated from fc-inflation model [1, H| . 
It is possible to regard the tachyon dark energy model [Iol - fl3T | , ghost condensation model [HI, EH , and scalar field 
quintessence model fl6rlT9| as variations of the fc-essence model. 

Since the fc-essence model is originated from fc-inflation model, it might be natural to consider a model unifying the 
inflation and the late acceleration by a single scalar field. In this paper, we try to construct such models by using the 
formulation of the reconstruction [20-25] and we also propose a mechanism of the reheating by the quantum effects. 
In the models, the solution which describes the unification of the inflation and the late acceleration can be stable in 
the space of solutions and also there does not appear tachyon due to the higher derivative terms. This tells that we 
can construct a model describing the time development, which cannot be realized by models of usual canonical scalar 
field like inflaton or quintessence due to the instability or the existence of tachyon. 

II. REVIEW OF THE RECONSTRUCTION AND THE STABILITY OF THE SOLUTION 

In this section, based on (25|, we review on the reconstruction by using e- folding N, which will be defined in this 
section, and discuss the stability of the solution in the space of solutions. A formulation of the reconstruction using 
the cosmological time has been given in [24j (about the reconstruction of the canonical/phantom scalar field, see 
[2fjl [27| and about the general formalism of the reconstruction, see [20M25| ) . In the formulation using the cosmological 



time 



_24j . it is troublesome and difficult to discuss about the stability of the solution when matters are included. In 
the formulation using the e-folding N, as long as the iV-dependence of the matters are known, which is often true as 
we will see, it is easy to construct a model where the solution is stable. 
We now consider a rather general model, whose action is given by 

S = [ <£xj—g (A ~ K (0, X) + L mattcl ) , X = d' 1 ^ . (1) 
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Here </> is a scalar field. Now the Einstein equation has the following form: 

^2 - l^Rj = -K (0, X) 9iiv + 2K X [cj>, X) d^ducj) + T M „ . (2) 

Here Kx {<t>,X) = dK (4>,X) jdX and T^ v is the energy-momentum tensor of the matters. On the other hand, the 
variation of <f> gives 

Q = -K^ > ( ( f>,X)+2^(K x (4>,X)d fi 4>) . (3) 



Here K$ (cp, X) = dK (<jf>, X) jd<f> and we have assumed that the scalar field <fr does not directly couple with the matter. 
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We now assume the FRW universe whose spacial part is flat: 

ds 2 = -dt 2 + a{t) 2 ( dxl ) 2 > ( 4 ) 

j=l,2,3 

and the scalar held cj> only depends on time. Then the FRW equations are given by 

Z_ H 2 = 2X dK< ^ x X) - K (0, X) + /Wter(t) , ~ (2H + 3ff 2 ) = K (0, X) + JWterW • (5) 

It is often convenient to use redshift z instead of cosmological time t since the redshift has direct relation with 
observations (see [23[ for the reconstruction of F(R) gravity using the redshift z). The redshift is defined by 

°W = 07) = eN ~ N ° ■ ^ 

Here to is the cosmological time of the present universe, No could be an arbitrary constant, and N is called as e-folding 
and directly related with the redshift z. In terms of N, the FRW equations (|SJ) can be rewritten as 

3_ h2 = 2X dK(^,X) _ R ^ + pmattoi(7V) ^ _J_ ( 2HH , + iH 2^ = K ^ x) + p mattcr (7v) . ( 7 ) 
Here H' = dH/dN. If the matters have constant EoS parameters Wi, the energy density of the matters is given by 

i i 

p mattei (N) = £ Wi p 0i a- S ^ = £ u^e" 3 ' 14 ^-^ . (8) 

Here poi's are constants. Eq. flSJ) tells the N dependence of the matter energy density p ma ttcr is explicitly given. Note 
that the ./V dependence is not so clear when the matters are created or annihilated as in the period of the reheating 
but in the periods of the inflation and the late acceleration, the expression of /o ma tter in (JSJ could be valid. For the 
general energy density of matters /3 ma tter(-^0i since the conservation law 

^matter + 3-ff (/^matter + Pmattcr) = , (9) 

can be rewritten in terms of N as 

PmattcrW + 3 (Pmatter(JV) + JWter(W)) = , (10) 

we find 

Pmattcr(iV) = - p matteI {N) ~ jitter (#) ■ (H) 

Then we can rewrite the FRW equations (J7J as 

1 / dH ,\ , , 1 , , , dK(6,X) 1 dH 1 . , , , , 

K {^X) = -— \2H— + 3H 2 j + p matteI (N) + -/4 attcr (A0 , -X ^ ' = -ff— - -p' mattei (N) . (12) 

If we define a new variable G(N) = H(N) 2 , the equations in ([T2j) have the following forms: 

K(cj>,X)=~ (G'(N) + 3G(N)) + /WtcrW + fatter W , = A G ' ( N ) ~ ^matter W . (13) 



K 



grmattcrv ^ 2r2 - v / gf 



Then by using the appropriate function g<p(4>), if we choose 



X(0,X) = V( J +1) ^WW), 

t^O \H{$) + Tftn»"erW / 

X<°)(0) = (<£(0) + 3<fo(0)) , #«(</>) = 1 g' U) , (14) 
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we find the following solution for the FRW equations ([5]), 

G(N) = H{Nf = g<j> (N) + yp roatte r(iV) , <fi = N (X = -H 2 ) . (15) 

Now KW(<I>) with n > 2 can be arbitrary. As we will see, K^ 2 \<p) is related with the stability of the solution and the 
existence of tachyon although K^((/)) with n > 2 does not affect the development of the expansion of the universe. 

We should note that the solution (fT5"]) is merely one of solutions of the FRW equations (JT2J) in the model given by 
(JH) . In order that the solution (fT5|) could be surely realized, the solution (fT5|) should be stable under the perturbation 
in the space of solutions of the FRW equations. We now write the perturbation from the solution f|15|) as follows, 



G(N) = Go(N) + SG(N) \G (N) = g<j>(N) + —p m a.ttet(N) J , cf> = N + 5</>(N) . (16) 

We should note that in many cases, the iV-dependence in the energy density p ma ttcr of matter is usually given by a 
fixed function as in ([5]) and therefore we find 5p m atter = 0. Then the equations in (|13p gives, 



Then we find 



5(j>'(N) \ _ 1 ( A B\( 5(j)(N) 
6G'(N) ) ~ L(N) \ C D ){ 5G(N) 



Here 



2G (N) v " ~ 2G (JV)' 
C*^ (. 9 ;'(iV)+3 g ;(iV))L(iV)+3 ff ;(7V) 2 , D = -3i(JV)-3^(JV). (18) 

L(N)^g^(N)-8K 2 K^(N). (19) 



In order for the solution (| 1 5[) to be stable, the perturbations Sc/>(N) and 8G(N) should decrease with the increase of 
N, which requires that the real parts of the eigenvalues for the matrix f ^ ^ ^ should be negative. Therefore the 
stability of the solution requires A + D < and AD — BC > 0, which gives 

(2G (N) + 3L(N)) 4{N) + (g'^N) - G' (N) + L(N)) g'^N) - L(N)G' (N) > . (21) 

We can find H < 3H 2 from (|2T)|) . which is always satisfied when the universe is in the non-phantom phase, where 
H < 0. For later convenience, we rewrite (f2~Tj) in the following form: 



H 2 K^'(N) - 2HH'K {1 \N) + 3K 2 k^\N)k^'{N) + 2 [K^(N) 
+ (iHH' - l2n 2 k {1) '{N) - 4 K 2 A" (1) (iV)) K {2) {N) > . (22) 

The condition (f2"Tj) or (|2"21 can be satisfied by choosing L(N) and therefore k^ 2 \N) properly. 

In case of usual inflaton or quintessence model, where k^ n \<p) — (n > 2) in (fT4|) . there appears tachyon if the 
potential is concave downwards and therefore the system becomes unstable. We now show that the development of 
the expansion in the universe generated by the concave potential in case of the inflaton or quintessence model can be 
realized without tachyon by adjusting k^ 2 \4>) in the fc-essence models in this paper. 
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We now consider the perturbation of only scalar field <fi from the solution (1151) as 

(j) = N + 8<j){x l ) . (23) 



Different from the case of (|16[) . we assume 8<f> only depends on the spacial coordinate x l since we are now interested 
in the pole of the scalar field propagator for the spacial momentum, corresponding to tachyon. Then by using §3§ and 
(fl4l) . we obtain 

n _ o* (1 W A ,^ , Ji^oVV.n , *<lY>r*n , (o H '\ rXlYr^ , ( , f H'\ 2 W> 



H 2 a 2 



A (6cP) + 2 W(iV) + KM"(N) + (3 - f ) ^ (1) '(iV) + I + - ^ j *T«(iV) 



iT\ 2 4ff" YlW \ 4H 7 - 



+ (^-4^J +^- + ^j^ 2 )(^) + ^ 2 )'(^)|^. (24) 
Here A is the Laplacian for the spacial coordinates x*. Then if 

L*kwm + (-- + f-V - — ) mm 

KW(N) \ H \ H \ H / H ) 

+ f-4 (Tf + %- + ^ ^ (2) W + ^> W < , (25) 



/ # H y ' H 

there does not appear tachyon. If we assume K^ X \N) < 0, which corresponds to non-phantom universe, and define 
U 2 \N) by 

(N) = H (^H'y 1 (N) , (26) 

Eq. ([25]) gives 



dU 2 



<p=N 



a 3 iH'~ nY , x / H" fH'V 6H'\ - m , A , . 

> T ^ (N) - -- + (-) - — . (27) 



Then if (I2lj) or (|22j) and (1251) or (|27[) are satisfied simultaneously without divergence, we obtain a stable model without 
tachyon. 

III. MODELS UNIFYING THE INFLATION AND THE ACCELERATING EXPANSION 

In this section, we propose models unifying the inflation in the early universe and the accelerating expansion in the 
present universe. 

In (fT5|) . g<f,(N) corresponds to the energy density p<j> of the scalar field 0: 

MAO = -| g0 (AO = 2X dK gx X) -K(4>,X) . (28) 

We expect that the energy density p<§, would behave as the cosmological constant in the period of the inflation and the 
late acceleration. Then we expect the behavior of p^ as in FIG. [T] We consider the model that the particle production 
and the reheating would occur after the inflation. 
We now assume 

1. The energy scale of inflation should be almost equal to the GUT scale. 

2. Except the period of the particle production, the EoS parameter w$ for the scalar field cj> could be given by 

-'-^m- (29) 



5 




FIG. 1: The expected behavior of p$(N). 



3. In general, there are two solutions N = AT l5 N 2 (N% < N 2 ) in the equation 

^ N=N 1<2 1 



(30) 

We expect that the expression ([29]) could become invalid when N\ < N < N 2 . 

4. The inflation started at N — and the end of the inflation is defined by N = N\ = Ni ~ 60. 

5. The reheating and the particle production could have occurred when N\ < N < N 2 . 

6. The reheating temperature T RH could be lOMeV < T RH < 10 14 GeV. 

Furthermore, the cosmological observations tell, at present, 1) the energy density of the dark energy is about 
10~ 47 GeV 4 , 2) the temperature of the present universe is 2.725 K and that at the epoch of the decoupling is al- 
most 3000 K (0.26 eV), which give the following constraints 

1. P4,(N = 0) ~ 10 60 GeV 4 , 

2. p (iV = jV o )~lO- 47 GeV 4 , 

3. w<j) = -1.023 ± 0.144 at TV (31) 

Here we choose No as the e-folding at present universe and the third constraint in (|3"Tj) comes from SuperNova Legacy 
Survey (SNLS) date [|. 

We should note that in the period of the reheating and/or particle production, it is difficult to apply the formulation 
of the reconstruction since the matter energy density is not always given by an explicit function of the e-folding N. 
In this paper, we approximate the behavior of p^, in the period of the reheating and/or the particle production by the 
interpolation from the behaviors in the period of the inflation and that after the reheating. 



A. Model 1 



We now consider the following model as model 1: 



P m = M 4 exp ( -— ri ^ x | , (32) 

d 1 + cicxp ( 
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which gives 



p'<f,(N) 



exp 



( n-nA 



ciA 



^( M )-i+exp(-^)) 



(33) 



Here c%, d, and Ai are constants and we choose c\ ~ 6.309 and M ~ 10 15 GeV. Then the assumptions mentioned 
above and the constraints (pn) give 




1 - 107 In 10 • ci exp 



Np-N! 
A! 



(iVo - JVi) 
Ai 

exp 



Ai 



0.363ciAi 

Since the scale factor a is proportional to the inverse temperature a = e N ~ N o cx T -1 , we find 

N _ Ni _ a(N) a(N) T RH 



a(Ni) a(N } 



RH J 



and therefore 



In 



3 x 10" 4 eV 



= 24-61 for T RH = 10 MeV-10 14 GeV . 



(34) 



(35) 



(36) 



The second constraint in (|34|) tells that the parameter d is expressed by the another parameter Ai, so in this model 
there remains only one undetermined parameter. Then in case Trh = lOMeV, we find (0,246.4) < (Ai,rf) < 
(1.81,247) and in case Trh = 10 14 GeV, (0,246.4) < (A u d) < (4.97,248.2). 
Now the reconstructed action has the following form: 



K(<f>,X) 



K {1) (<f>) = -AT 4 cxp - 



Eh- 




X 



3 P<t>(<t>) + ! TPm(4 l ) 

exp 



+ 1 K {n) (cj>), 



N-Nj 



\ 



3ciAi ((cid)- 1 + exp(- 



N-Ni 



J 



exp 



2 • 



In order to find the constraints for K^ 2 \<p) or k^ 2 \<p) given by (|21|l or 

p m (N) 



edAifM^+exp^-^p 

and (|2"5j) or (|2"T1) . we assume 



(37) 



PmO 



I for < TV < ATi 

~ \p m oe- 4{N - No) for JV > ATrh * N 2 , 
8.4 x 10" 52 GeV 4 . 



(38) 



Here -/Vrh expresses the e-folding number when the reheating finished. Then, we obtain approximate constraints for 
k^ 2 \(j)) and k^'(cj)) for model 1 as shown in FIGs. HH7] and we can find that there exists k^ 2 \(f)) or K^((/)) which 
satisfies the constraints and does not have divergence nor vanish. 

In FIG. [H the region satisfying the constraint (|2"2")l is depicted by the directions of arrows. When N\ = N% < N < 
N2 ~ -/Vrh, we could not be able to use the formulation of the reconstruction due to the particle creation. The region 
< N < Ni in FIG. [2] is magnified in FIG. [3] and the region N > N 2 ~ N RH in FIG. H The region that k^'(N) of 
model 1 satisfies the constraint ([27]) is depicted in FIG. 0] and regions < N < Ni and N > N 2 ~ jVrh in FIG. 0] 
are magnified in FIG. [5] and FIG. [51 respectively. Then we may find that we can always obtain an action where the 
solution becomes stable and does not have tachyon. 
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FIG. 2: The regions satisfying the constraint ()22[) for k^ 2 '(N) of model 1. The gray regions express the forbidden regions for 
the instability of the solution. In the interval Nx = Ni < N < N2 — Nnu, the formulation of the reconstruction could not be 
applied due to the particle creation. 




FIG. 3: The region < N < Ni in FIG. [2] is magnified. The vertical axis expresses the absolute value of k^ 2 '(N). The symbol 
'n' means the value of k (2) (N) is negative there. 




70 80 90 100 110 



FIG. 4: The regions N > ATrh in FIG. are magnified. The vertical axis expresses the absolute value of k^(N). The symbol 
'p' means the value of k^ 2 '(N) is positive there. 
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1.0x10' 




100 120 



FIG. 5: The regions satisfying the constraint (|27l) for k^'(N) of model 1. The gray regions express the regions forbidden by 
the constraint (|27[) . 





FIG. 7: The regions N > Nbm in FIG. [5] are magnified. The vertical axis expresses the absolute value of k^ 2 ''(N). 
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B. Model 2 



As a second model, which we call as model 2, we consider the following: 



P<t> (N) = j^— +B{N + b)- p , (39) 

c 2 exp(^)+l 



which gives 



Ac 2 exp(^f) 



p'^N) = — ^ " 2 ' - -f—B (N + by . (40) 



(c 2 exp(^) + 



N + b 



Here c 2 , A, B, A 2) and b are constants and we choose c 2 ~ 0.114 and A ~ 10 60 GeV . We now assume that the term 
with a coefficient A would dominate in the expression of p$ in (|3"Tj) when < N < Nj and the term with the coefficient 
B would dominate when N > _/V 2 . Furthermore we choose A — B in order to reduce the number of parameters. Then 
the constraints (|3"Tj) give 

2. 

107 In 10 



3. /? 



ln(JV + b) 



The second constraint in (|4ip tells that the parameter 6 is expressed by another parameter /3, so this model has two 
undetermined parameters /3 and A 2 . Then in case T RH = lOMeV, we find (0,0) < (A 2)/ S) < (0.095,47.26) (6 > 99.6) 
and in case T RH = 10 14 GeV, (0, 0) < (A 2 , /3) < (0.241, 49.01) (6 > 31.5). 
Now the reconstructed action has the following form: 

k ^ = 2 ( !T+ <^ x+k{a) + k(i)+ fi{ i? ,^ x * 7^ +i ) ^ (n) W' 



'■2 



* W W = " I 1 " ' I ^ ( 1 - 8^ ^ ( ^ + ^ 




KWtf) = °- V 2 \i -7rrT-zM<l> + b)- P ■ (42) 



(c 2 exp(^)+l) 



6(0 + 6) 



Similar to the model 1, by adjusting K^{(j)) of model 2 which does not have divergence nor vanish, we obtain a stable 
model without tachyon. In FIG. [51 the region satisfying the constraint (|2"2"]l is depicted and the region satisfying the 
constraint ((27)) is depicted in FIG. |H 



C. The dynamics of the scalar field 

The evolution of the expansion in universe does not change even if we consider the model with = (n > 2), 
which corresponds to the usual inflaton and/or quintessence models since the time evolution of the system is controlled 
only by K^((f>) and {</>). In case of K( n > — (n > 2), the scalar field becomes canonical and the dynamics of 
the scalar field is compared with the dynamics of a classical particle in a potential. In case of KM = ( n > 2 ), m 
order to generate the development of the universe expansion given by in the previous Subsections IIII Al and IIII B[ the 
potential has typically the form depicted in FIG. [TD] 

As an initial condition, the scalar field should almost stay near the top of the potential in order to generates the 
inflation. After that, it rolls down to the bottom of the potential and creates the particles. Finally, without getting 



Is 
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FIG. 8: The regions satisfying the constraint ([22]) for fe (2) (N) of model 2. 
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FIG. 9: The regions satisfying the constraint (J27J for k (2) '(N) of model 2. 




FIG. 10: The effective potential of the canonical scalar field. 
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trapped in the bottom of the potential, the scalar field goes through the subsequent small peak of the potential and 
plays the roll of the dark energy. 

It is important that different from the inflaton or quintessence models, we need not to fine-tune the initial conditions 
for the scalar field and there are no tachyonic instability in the models we have constructed even if the effective potential 
is concave downwards since the motion of the scalar field can be stabilized by their K^ 2 \4>) term which should not 
vanish. 



IV. A MECHANISM OF THE PARTICLE PRODUCTION 



Now we assume the Hubble rate is given in terms of the e- folding N as H = H(N) and consider the situation that 
the e-folding N can be identified with a scalar field <fi. We now consider the interaction between the scalar field <p 
between another real scalar field ip as follows, 

C ° f j3 / — d M</0 2 



ffi »--Y/ d ^T^- (43) 

Here Co is a constant. Note that p (/>((/>) is not the real energy density of </> but merely a function of <fi given by replacing 
N in p<p(N) in (12511 by cp. We assume that cj> can be treated as an external source and the interaction occurs only 

in a narrow region around t = and we approximate Cq as a function of the cosmological time t. We now 

approximate Cq dp %jjf* by the Gauss function: 



(44) 

Here Uq is a constant and A is the standard deviation. We also assume that the space-time can be regarded as static 
and also fiat when \t\ ~ A, which should be checked. 

Then the amplitude that the vacuum could transit to two-particle state whose momenta are given by p and q is 
given by 



/oo 
dt(p,q\H- mt \0) 
-oo 



TT fOO n -^-i(u p +uj q )+i(p+q)-x 

" 1 dt d 3 x- 



l S 3 (p + q)^- . (45) 



Here ui p = y p 2 + m^, with the mass m v of ip. Then the transition probability is given by 

1 f ,3 ,3 c3 / ,2 U 2 e~ 2A2 < VU 2 f 2 , e~ 2A2 < 



The factor 1/2 in the first line appears since (p,q\ = (q,p\ and V is the volume of space which appears since 

V_ 

(27) 

Especially when tp is massless, that is, m v — 0, we find 



s3 0) = — ■ («) 



8(2tt) 5 A 



which diverges when A — > 0. 

Then the total transition probability per unit volume is given by 



P2 = y- ( 48 ) 
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Therefore the particle density n is given by 

n = 2 P2 . (49) 

We now consider about the energy (density) . Eq. (J45J tells that the expectation value of the energy E 2 corresponding 
to two particles state is given by 

1 f *3 * / , ^ 2^C/ 2 e- 2A2 ^ VU 2 n [ , , e~ 2A2 



E 2 = - j d*pd* q 5* (p + g)' = ^ / ,r,/y, 



Especially when (p is massless, we find 



E 2 = -^f- . (50) 
8(2tt) 5 A 

Then the expectation value of the energy density e 2 for the two particle state is given by 

e 2 = f- (51) 
We may estimate the width A in ([HJ) by using Ni and ./V2 in (1501) as 

AiV = iV 2 - iVi = I Hdt~Hi f dt= 2#iA , (52) 



which gives AiV ~ 2.98Ai for model 1 and AN ~ 4.34A2 for model 2. Then since the energy density of the radiation 
in the present universe is given by the product of the critical density p cr o and the density parameter Sl r o for the 
radiation. Since 

Pcro = 1(T 47 GeV 4 , n r0 = 8.4 x 10~ 5 , (53) 

we find 

£2 . . 10.4 x 10-»-10«GeV - 4 ;~^ , (54) 

which tells 

U = 9.86Ai x 1(T 15 -10 17 GeV for model 1 , 1.43A 2 x 10" 14 -10 18 GeV for model 2 . (55) 
In (I54p . Tq is the temperature of the present universe. 



V. SUMMARY 



In this paper, after reviewing the formulation of reconstruction for fc-essence, we explicitly constructed two models 
which unify the inflation in the early universe and the late-time acceleration in the present universe, and satisfy the 
observational constraints. We have proposed a mechanism of the interaction for particle production by the quantum 
effects coming from the variation of the energy density of the scalar field and estimated the energy density of the 
particles. 

In both of the models, the solutions describing the development of the universe expansion are stabilized by or 
terms which should not vanish. We also note that K^ 1 or k^> terms play the role to eliminate the tachyon. As 
explained in Subsection IIII Cl the solutions describing the development of the expansion in our models can be realized 
by the usual inflaton or quintessence model, where the scalar field is canonical, but in the canonical scalar models, 
the solutions could be often unstable and there could appear a tachyon when the scalar field lies at the concave part 
of the potential. The instability of the canonical scalar models require the fine tuning of the initial conditions, which 
makes the models unnatural. In our models, due to the stability of the solutions controlled by the or fc' 2 ' terms, 
there could exist a wide region of the possible initial conditions. Then in the framework given in this paper, we 
can construct a model describing the time development, which cannot be realized by a usual inflaton or quintessence 
model. 

The roles of (n > 3) are, however, still unclear although these terms play the role to guarantee the existence 
of the Schwarzschild solution [22, 25]. More detailed cosmological constraints may restrict the form of these terms. It 
might be interesting to consider the reconstruction of the general spherical symmetric solution in the fc-essence model 
as in [H. 
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